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The role of Lorentz invariance as a fundamental symmetry of nature has been lately reconsidered in 
different approaches to quantum gravity It is thus natural to study whether other puzzles of physics 
may be solved within these proposals. This may be the case for the cosmological constant problem. 
Indeed, it has been shown that breaking Lorentz invariance provides Lagrangians that can drive the 
current acceleration of the universe without experiencing large corrections from ultraviolet physics. 
In this work, we focus on the simplest model of this type, called GCDM, and study its cosmological 
implications in detail. At the background level, this model cannot be distinguished from ACDM. 
The differences appear at the level of perturbations. We show that in GCDM, the spectrum of CMB 
anisotropics and matter fluctuations may be affected by a rescaling of the gravitational constant in 
the Poisson equation, by the presence of extra contributions to the anisotropic stress, and finally 
by the existence of extra clustering degrees of freedom. To explore these modifications accurately, 
we modify the Boltzmann code CLASS. We then use the parameter inference code Monte Python 
to confront GCDM with data from WMAP-7, SPT and WiggleZ. We obtain strong bounds on 
the parameters accounting for deviations from ACDM. In particular, we find that the discrepancy 
between the gravitational constants appearing in the Poisson and Friedmann equations is constrained 
at the level 1.8%. 



I. INTRODUCTION 

Explaining the origin of the current acceleration of the 
universe is one of the biggest challenges in cosmology. 
Despite the great successes of the ACDM paradigm, it is 
hard to accept that a very fine-tuned form of otherwise 
invisible energy governs nowadays the behavior of the 
universe at the largest scales. A whole plethora of mod- 
els are proposed as alternatives to this situation. These 
are known as quintessence or dark energy models [IH3]. 
The motivations and domain of applicability of various 
proposals are quite diverse: whereas some of them are 
intended as mere phenomenological models only valid for 
cosmology, others are rooted in theoretical considerations 
and are testable by different types of experiments. From 
the theoretical viewpoint the preferred models are those 
addressing (at least partially) the naturalness problem of 
the cosmological constant without introducing any addi- 
tional fine-tunings. In this paper we will study in detail 
the OCDM proposal of Ref. [4] that has such a property. 
This model is well-motivated theoretically and has a rich 
phenomenology that may distinguish it from ACDM. 

The OCDM model is based on the idea that the 
Lorentz invariance observed in the Standard Model of 
particle physics is an emergent phenomenon and does 
not correspond to a symmetry of nature. This concept 
is motivated by attempts to find complete theories of 
quantum gravity, but it can also be considered indepen- 
dently. There exist two different approaches: on one 



hand Einstein-aether theory [5j|6] provides a phenomeno- 
logical description of gravity with broken Lorentz sym- 
metry at large distances; on the other hand Horava grav- 
ity [7J [8] invokes a violation of Lorentz invariance at any 
scale to improve the quantum properties of gravitational 
theories. The restriction of Hofava gravity to operators 
with the lowest dimension can be considered as an effec- 
tive theory by itself, called the khronometric theory [9]. 
The two proposals are very similar at large distances: the 
khronometric theory can be viewed as a constrained ver- 
sion of the Einstein-aether theory, and in many situations 
the predictions of the two theories coincide [9j [10] . The 
analysis presented in this work is applicable to all these 
classes of models. However, in the last step consisting in 
deriving the actual observational bounds on the model 
parameters, we will restrict to the khronometric case. 

Several implications of the above models for cosmology 
have been studied in the past, see e.g. [TTlfT7] . For the 
evolution of the Universe as a whole, the consequences of 
the minimal models 1 are rather trivial. They are compat- 
ible with Friedmann-Robertson- Walker (FRW) solutions 
that differ from the standard case only by a renormaliza- 
tion of the gravitational constant away from the value 



By this term we refer to the models with a minimal number of 
additional degrees of freedom in the gravity sector, and the sim- 
plest kinetic action (Einstein-aether and khronometric theories) 

E][9]. 
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measured in local tests of Newton's law. In other words, 
those models do not modify the form of the Friedmann 
equation. The background evolution becomes more in- 
teresting when non-minimal models are considered. In 
this case, the breaking of Lorentz invariance allows for 
Lagrangians which can change the expansion history of 
the universe and provide new alternatives for inflationary 
dynamics, dark matter and dark energy [4j [TTJ [T71 [18]. 
Some criteria are necessary to identify the most interest- 
ing cases. In this work we will be concerned with the 
issue of dark energy, for which one would like to find a 
simple model described by a Lagrangian with high cutoff 
scale, that provides a mechanism to accelerate the uni- 
verse insensitive to UV corrections 2 and distinguishable 
from ACDM. The 6CDM model of Ref. @] meets these 
requirements. 

Some consequences of BCDM have been already dis- 
cussed in Ref. [4 . In the present work we describe the 
observable physical effects of the model on cosmic mi- 
crowave background (CMB) anisotropics and on the mat- 
ter power spectrum (at the linear level). We provide a de- 
tailed discussion of such effects, that we computed accu- 
rately with a modified version of the flexible Boltzmann 
code CLASS [19]. We then compare the 6CDM model to 
recent CMB and Large Scale Structure (LSS) data us- 
ing the Monte Carlo parameter inference code Monte 
Python 20]. 

Our work is organized as follows: in Sec.[ll]we describe 
the OCDM model and discuss the constraints not re- 



lated to cosmology. In Sec. [Hi] we study the background 
evolution and derive linear equations for perturbations 
around the FRW background. The results presented in 
this section are complementary to those in Ref. [4 . They 
hold in a different gauge, and refer to different conven- 
tions and parametrizations, found to be more suitable for 
the numerical implementation. The qualitative effects of 
OCDM on the CMB and matter power spectrum are de- 
scribed in Sec. |IV| We present constraints from CMB and 
LSS data in Sec.|V[ and expose our conclusions in Sec. |VI| 
Appendix [A] contains the derivation of initial conditions 
for cosmological perturbations in OCDM. 



ity to a complete quantum theory, which points towards 
the necessity to consider alternatives. Independently of 
this, modifications to general relativity are currently be- 
ing considered in the area of cosmology. The rationale 
behind these modifications is the possibility to use the 
wealth of cosmological data to learn how gravitation be- 
haves at the largest accessible distances and, hopefully, 
shed some light on the mechanism responsible for the 
accelerated expansion of the universe. 

These two lines of research converge if one assumes 
that Lorentz invariance is not a symmetry of the gravita- 
tional sector. This idea opens the possibility to construct 
gravitational theories with better quantum behavior than 
general relativity [7J. This is achieved at the price of in- 
troducing new degrees of freedom that modify the laws 
of gravity at all distances, including those relevant for 
cosmology [8j [9]. Even before this top-down approach 
had been initiated, the bottom-up Einstein-aether model 
[5] was proposed as a way to capture the large-distance 
effects of a putative Lorentz violation due to quantum 
gravity. In this work, we will consider theories where 
Lorentz violation is described by a preferred time-like 
vector field defined at every point of space-time, which 
includes Einstein-aether theory and Horava gravity. For 
the latter, we will focus only on its low-energy form, the 
khronometric theory [9] . The vector field will be normal- 
ized to 3 



UnW 



-1. 



(1) 



The presence of this dynamical field allows to use for the 
description of Lorentz violation the same language as for 
spontaneous symmetry breaking. Gravitational physics 
at large distances is governed by the following covariant 
action for g^ v and featuring a minimal number of 
derivatives: 



M 2 (' 



(2) 



where 



II. LORENTZ BREAKING THEORIES OF 
GRAVITY AND OCDM 

Lorentz invariance is one of the best tested symmetries 
of the Standard model of particle physics [21] . It is also 
a fundamental ingredient of the theory of general relativ- 
ity, that provides a very successful description of gravita- 
tional interactions over a huge range of scales. However, 
it is not known how to directly promote general relativ- 



K»\ p = c ig ^g ap + c 2 5^ p + c^S" - c 4 u^g ap , (3) 

and / is a Lagrange multiplier that enforces the unit-norm 
constraint. Equation ([2| includes the Einstein-Hilbert 
(EH) term evaluated with the metric g^ v . The parameter 
Mo is proportional to the Planck mass (cf. (|8|) while the 
dimensionless constants c a , a = 1, 2, 3, 4, characterize the 
strength of the interaction of the aether with gravity. 
This is the action of the Einstein-aether model [5] [6] . 



2 We leave aside the "old cosmological constant problem": to find 3 We use the ( — h ++) signature for the metric. This differs from 
a mechanism imposing a null vacuum energy (see however the most of the previous works in the field, but is common in cos- 

related comments in [4].) mology. 
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The khronometric case corresponds to the situation 
where is defined as a vector normal to a foliation con- 
sisting of the level surfaces of the khronon field (/?, 



(4) 



In this case, the constraint is satisfied identically and 
the first term of ([3| can be expressed as a linear combi- 
nation of the last two terms. Thus, K^ v ' ap reduces to its 
last three terms with coefficients 



A = c 2 , /3 = c 3 + ci, 



c 4 + Ci . 



(5) 



At the level of linear perturbations, the khronometric 
theory differs from the Einstein-aether theory by the 
number of propagating degrees of freedom apart from 
the spin-two mode of the graviton. While the aether 
in general describes vector and scalar excitations, the 
khronometric case only contains scalars. However, the 
scalar sectors of the two theories are equivalent and are 
completely characterized by the three parameters 

Once the couplings to matter are specified, the con- 
stants c a are constrained by various considerations rang- 
ing from theoretical requirements to observational tests. 
First, to very good precision, the Standard Model fields 
must couple only to g^ u and not to , as required by 
Lorentz invariance in this sector 4 [21 . This decoupling 
presents a serious challenge to the proposal, but it is 
conceivable to achieve it either by imposing extra sym- 
metries, e.g. supersymmetry [23[[24], or through a renor- 
malization group running [25 j. Next, there are restric- 
tions imposed by the stability of Minkowski spacetime 
[6]. In particular, the requirement that the scalar mode 
is neither a ghost nor a tachyon field amounts to the con- 
straints 



0< a < 2 , /3 + A > . 



(6) 



Stringent bounds come from observations of the Solar 
System dynamics, which can be used to place constraints 
on post-Newtonian (PPN) parameters. Two of these pa- 
rameters, denoted by a p PN and a PPN , describe the ef- 
fects of Lorentz violation. In general, they are different 
from zero (their value in general relativity) both in the 
Einstein-aether and khronometric theory; we refer the 
reader to [6j [29] for explicit relations with the constants 
c a . Assuming no cancellations in these formulae, one ob- 
tains experimental bounds of the order of 



< 10" 



(7) 



However, there are regions in the space of model param- 
eters (a, /3, A) in which ol ppn vanish and Solar System 



4 More generally, a universal coupling of Standard Model fields to 
a fixed combination of and is allowed. This reduces to 
the previous case by a redefinition of the metric |22| . 



tests are automatically satisfied, with Newton's constant 
given by 



G 



N 



8ttM 2 (1 -a/2) 



(8) 



This requires 5 



a = 2fi for the khronometric model, (9a) 
a = — (3A + j3) for the Einstein- aether. (9b) 

In these cases, much weaker bounds can be inferred from 
gravitational wave emission in binary systems [28l [29] . 
giving 



< 10" 



(10) 



and from the form of black hole solutions, imposing in- 
equalities detailed in [26]. We will see that in the case of 
the khronometric model, the bounds that can be inferred 
from cosmology are competitive with ([l0| ), but not with 
0. Thus we will impose the relation rt9al) when searching 
for the allowed parameter space. Previous studies of the 
cosmological effects of Lorentz violation [TTJ [16] focused 
on the Einstein-aether model. In those works, the rela- 



tion (9b) was enforced to avoid the PPN constraints. As 



we will explain, the relation (9b) also incidentally sup- 



presses the leading effects in cosmology, unlike the re- 
lation ( |9a| ). This explains why the bounds obtained in 
those studies are rather mild. 

The implications of the action ([2| for the background 
evolution of an homogeneous and isotropic universe are 
minimal. If all matter components are universally cou- 
pled 6 to g^v, the only difference with respect to gen- 
eral relativity is that the Friedmann equation involves a 
renormalized gravitational constant 



1 



8ttM 2 (1 + /3/2 + 3A/2) 



(11) 



differing from the value Gn measured e.g. on earth or in 
the Solar System 7 , given by Eq. In order to modify 
the expansion history of the universe, and find a candi- 
date for dark energy, we need to add a new ingredient to 
the model. We want to do it in a way that preserves the 



5 The difference between the models is due to the contributions to 
a i2 N fr° m the vector polarizations that are present only in the 
Einstein-aether case. 

6 A priori, there is no reason why this must be true for the dark 
matter. Still, as shown in |27| , even allowing for non-universal 
interactions between the dark matter and the aether does not 
change the background evolution. As our main focus in this pa- 
per is dark energy, we will assume that the dark matter has stan- 
dard properties (namely, that it is a pressureless fluid, described 
at the fundamental level by a Lorentz invariant Lagrangian, and 
with universal coupling to g^v)- 

7 The analysis of Big Bang Nucleosynthesis 22 se ^s a bound on 
the relative difference between the two, \G CO s/Gn — 1| < 0.13. 
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simplicity of the proposal and its validity as a low-energy 
effective field theory. This is achieved by supplementing 
the action ([2| with a new field invariant under the shift 
symmetry, 



6^6 + const. 
The low-energy action for this field is 8 



(12) 



Si 




(13) 



and involves two free parameters (k, /jl). We will refer to 
the model resulting from the combination of the actions 
(13) and ([2| (with a universal coupling between matter 
fields and g^ u ) as 6CDM [4]. The cosmological constant 
term is set to be exactly zero 9 . Two comments are in 
order. First, as an effective field theory, OCDM provides 
a valid description of physics up to a cutoff scale of the 
order of A c ~ ^/c^M P . The bounds Q or ^0) show 
that this cutoff scale may be only a few orders oi magni- 
tude below the Planck mass. Furthermore, in the khrono- 
metric case, the theory has a potential UV completion, 
since it is a sub-case of Hofava gravity [7]. Second, in 
the whole OCDM action, only the last operator in (13) 
breaks the discrete symmetry — >> —0. This implies 
that from the viewpoint of the effective field theory, it 
is self-consistent to choose the coefficient in front of this 
operator to be much smaller than the UV cutoff. In spite 
of \i being dimensionfull, the above symmetry guarantees 
that it is renormalized multiplicatively, and that no dan- 
gerous contributions proportional to A c appear. In other 
words, the smallness of \i is technically natural. This last 
observation is very important, since we are going to see 
that in the 0CDM model, \i sets the scale of the current 
cosmic acceleration. 



III. COSMOLOGICAL SOLUTIONS OF OCDM 

We are interested in describing the evolution of per- 
turbations around homogeneous and isotropic solutions 
in the 0CDM model. We focus on scalar perturbations, 
and refer the reader to jT3l [TSJ [30] for possible effects of 
vector and tensor modes. In the synchronous gauge, the 



8 A similar action wit h an additional potential term for breaking 
the shift symmetry |l2| was considered in [18) . 

9 One can entertain the possibility of finding a mechanism that 
would enforce the cancellation of the vacuum energy induced by 
quantum loops |4]. However, at present, we are not aware of 
any such mechanism. Thus, the vanishing of the cosmological 
constant should be taken merely as an assumption. 



perturbed FLRW metric has the form, 



#00 



-a(r) 2 , #oi=0, 



9ij = «0) 2 



(14) 



Besides, for and the khronon (or the longitudinal com- 
ponent of in the more general Einstein-aether case) 
we introduce 



= 0(r)+£, ^ = t + X - 



(15) 



The matter components are assumed to be cold dark mat- 
ter (cdra), photons (7), neutrinos (y) and baryons (b). 
We describe these components at the same level of ap- 
proximation as in Ref. [31]. In particular, the dark mat- 
ter is treated as a pressureless fluid universally coupled 
to the metric g^ v . We assume that it is comoving with 
the gauge, in order to eliminate the well-known resid- 
ual freedom in the synchronous gauge. We now discuss 
how the gravitational equations are modified in 0CDM. 
These equations were derived in [4] using the conformal 
Newtonian gauge. Here we will rewrite the linearized 
equations in the synchronous gauge ( (14] ), and in a form 
optimized for numerical study with the Boltzmann code 
CLASS [T9] , 



A. Background evolution 



Deriving the equation of motion for from (13), one 



finds that the homogeneous part Q(t) evolves as 







2 



C 
^2 



(16) 



where C is an integration constant. Substituting this 
solution into the Friedmann equation (derived from the 
combination of the actions ([2| and ( |l3| )) yields 



H 2 = 



8ttG c 



P» 



other 



(17) 



where H = a /a and G cos is given by (11). The first 



three contributions in the brackets come from the energy- 
momentum tensor of the ©-field and have the form, 



2(1 + k) 



C 2 (l + k) 
2a e 



Pd 



[i 2 C 



(18) 

while p n , n = cdm,7, ^, 6, stand for the densities of the 
standard matter components of the Universe. 



Let us analyze Eq. (17). Recall that there is no cosmo- 



logical constant at the fundamental level in the model. 
Instead, the first term in (17) plays the same role with an 
energy scale set 10 by \i. As emphasized above, /i does not 



10 Throughout the paper we assume that the combination 1 + k is 
of order one. 
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receive large radiative corrections, and thus this source 
of dark energy can naturally have an energy scale com- 
pletely unrelated to the cutoff of the theory. The second 
term has the form of the contribution of stiff matter. 
Not to spoil Big Bang Nucleosynthesis (BBN), p s must 
be smaller than about 30 times the density of radiation 
in the Universe at temperatures of the order of 10 MeV 
[32] . Thus, assuming that this contribution was already 
present at BBN 11 , due to its rapid decrease, it is com- 



pletely negligible at later epochs. The third term in (17) 
behaves as the energy density of dust 12 , so one might be 
tempted to identify it with the dark matter. However, 
being the geometric mean of the first two, this term is 
always subdominant and cannot contribute a significant 
fraction of dark matter. Given these considerations, we 
will set C = henceforth. 



B. Cosmological perturbations 

Let us introduce two time scales that appear in the 
analysis of the linear perturbations, 



87rG r 



a 



e, H a 



Ho 



(19) 



where Hq = 100 h km s _1 Mpc -1 is the current value 
of the Hubble constant today. It is also convenient to 
rescale the O-fluctuation defining 



where we introduced the notations 



G = 



B 



1 


G 'cos 


8ttM 2 ' 




/? + 3A 




1 

a 


4- 



1 + /3/2 + 3A/2 ' 



a 



4 



and 



H = aH. 



(22) 



(23) 



Note that the constants c x , cq represent the sound speeds 
of the fields x an d £ at relatively short wavelengths, 
where these modes decouple from each other [4]. These 
two velocities are constant in time and we will treat them 
as quantities of order one, which is their natural order of 
magnitude from the point of view of effective field the- 
ory. Note that in general they can exceed one, so that the 
fields can be superluminal. In Lorentz violating theories, 
this does not lead to any causal paradoxes, see e.g. the 
discussion in Ref. [33] . 

We now present the linearized Einstein equations. 
There are four equations in the scalar sector correspond- 
ing to different components of the Einstein tensor G^ v = 
W v — \5^ V R. Only two of these equations are indepen- 
dent, but we write all of them for completeness. For the 
SG° component we find 



^ ^COS 



-WGo PnSn, (24) 



with 



~ V87rG cos 

S — TT S 5 

^0 



(20) 



so that both fields £) have the dimension of time. The 
equations of motion for x an d £ read 



X = ~ 2Hx 



k 2 c 2 + (1 + B)U 2 + (1 - B)H 



Go 



G T 2 
^cos ' a _ 



Go Ha z C x I o ' • 



G C os 7~a 



P. 

—r 

a 



(21a) 
(21b) 



11 We do not specify the origin of the field in this work. One 
possibility is to identify with the Goldstone boson of a spon- 
taneously broken global symmetry. Then, the above assumption 
amounts to stating that the corresponding phase transition oc- 
curs before BBN. If it happens later, the picture may change, 
but it is hard to see how this option can be incorporated in a 
viable cosmological scenario. 

12 Note though, that its sign can be negative depending on the sign 
of C. 



= ^2 pA 



87ra 2 G cos ck r a 



other u ^5) 

Here and in what follows, the label 'other' refers to con- 
tributions from the standard matter components, whose 
form can be found in [31]. For the diSG° i part we find 



2k 2 (l - 0)ri - -^k 2 h = 8™ 2 G Y,(p n + Pn)0 n , (26) 

n 

with 



a c 



other 



87ra 2 G 



k\. (27) 



The SG l i equation reads 

.G c 



h = -2Uh + 2-^k 2 rj - 2^G cos a 2 V Sp n , (28) 



where 



aB 



other 



247ra 2 G 



k 2 (x + 2U X ). (29) 
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Finally, the didjSGj equation yields 



where 



(1 -/3)(h + 6rj + m(h + 6?))) - 2k 2 T] = 

-24na 2 G ^(pn + Pn)v n , (30) 



with 



/3k 2 



(Pn +Pn)CT n = V (p n +Pn)<J n - 

-(X + 2WX)- 



other 



127ra 2 G 



(31) 

In the above expressions the pressure fluctuations 5p n , 
the velocity divergences n and the shear potentials a n 
are defined in the standard way, see [31]; the equations 
for perturbations in the matter components can be found 
in the same reference. The linearized equations must be 
supplemented by suitable initial conditions. The latter 
are derived in Appendix [A] 

For simplicity, we are going to set cq = 1 in the nu- 
merical simulations. This prescription is equivalent to 
fixing k = 0, and leaves us with four free fundamental 
parameters (a, /?, A, /i). A different choice for cq would 
not affect qualitatively the evolution of perturbations, 
unless cq is very small. In that case, the field £ could in 
principle cluster on scales much smaller than the Hub- 
ble radius, but we will not consider this situation in the 
present work 13 . 

We can compare the number of free parameters in this 
model with that in an ordinary ACDM model sharing the 
same background evolution. The parameters (ft a, Hq) of 
ACDM can be mapped onto the parameters (/i, Ho) of 
the 0CDM model. Hence the latter features only three 
additional parameters (a, /?, A), reducing to only two 
independent parameters after imposing one of the condi- 
tions fulfilling Solar System tests, (9b) or (9a). Once (a, 
/3, A) are fixed, all coefficients in the field equations (21) 
and in the Einstein equations can be derived. 



C. Changing the gauge 

Reference [31 shows the impact of a gauge transforma- 
tion on the variables describing matter and metric pertur- 
bations, and on the form of their respective equations of 
evolution. Here we show how the fields x an d £ transform 
when switching from the synchronous to the Newtonian 
gauge. This transformation is induced by the particular 
change of coordinates 

x° i-> x° - $(x, r) , x* i-> x* - dif3(x, r) - e\x) , (32) 



P(x,t) 



I< 



dJ = . 



(33) 



After this transformation, the scalar part of the metric 
(14) becomes diagonal, 



ds 2 = a(r) 2 [-(1 + 2^)dr 2 + (1 - 2(j))dx i dx i } , 
with 

^ = ^ + H/3, (j) = rj-HP . 
The khronon and field transform as 



(34) 



(35) 



N 



Ho 



-6/3. (36) 



Although the Boltzmann code CLASS features equations 
in both the synchronous and Newtonian gauge, for sim- 
plicity, we choose to implement the 6CDM equations 
in the synchronous gauge only. However, when present- 
ing the physical interpretation of numerical results in the 
next section, we will refer to the evolution of quanti- 
ties in the Newtonian gauge, obtained by performing the 
above transformation inside the code a posteriori (i.e. 
after solving the equations of motion in the synchronous 
gauge). 



IV. OBSERVABLE EFFECTS 

The modified evolution of linear perturbations in 
OCDM leads to observable consequences which we now 
discuss. We focus on the scalar sector, known to 
contribute most to observable quantities. The effects 
of Lorentz violation on the tensor and vector sectors 
are weakly constrained by current cosmological data 
[T3l [30] . The only effect on tensor modes is a small 
shift in their velocity [6 j. The vector equations are iden- 
tical in general relativity and in the khronometric model, 
implying that vector perturbations decay with time. Pos- 
sible effects of vector modes in the Einstein-aether model 
on CMB polarization have been discussed in p~3j [15] . 

At the qualitative level, one can identify three main 
effects distinguishing the growth of perturbations in 
OCDM from that in ACDM. These are: (i) a rescaling of 
the matter contribution in the Poisson equation (i.e., a 
different self-gravity of matter perturbations), (ii) an ad- 
ditional contribution to the anisotropic stress, and (in) 
the presence of additional clustering species. The first 
two effects are generic for any Lorentz violating gravita- 
tional theory based on the Einstein-aether or khronomet- 
ric model 14 , while the third is specific to the dynamical 
realization of dark energy in OCDM. 



A small value of cq corresponds to the near cancellation of the 
kine tic t erms for £ coming from the first and second terms in 
Eq. ( |13| - a situation that appears fine-tuned from the effective 
field theory perspective. 



Similar effects are also known in other modified gravity and dark 
energy models, see e.g. 34 35 . 
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To understand these effects, let us work in the Newto- 
nian gauge. The Poisson equation (i.e., the sub-Hubble 
limit of the (00) Einstein equation) reads 



/c 2 (20-a?/O 



other 



a G 



ak 2 (x + m ~ B) X ) ~ 2 ( H a ti 



t 

T„. 



(37) 



Let us first concentrate on the contribution of standard 



matter. Using the Friedmann equation (17), the (time- 
dependent) fraction of the total energy density of the 
Universe due to each matter component is given by 



In — 3^"2 P n ' 

The Poisson equation takes the form 

k 2 cj) = -47rG N a 2 ^ p n S n 
3 Gn 



(38) 



other 



2G C 



-7~L 2 fn$n 



other 



where for the time being we have omitted the terms in the 
second line of (37) and neglected the anisotropic stress 
(i.e. we have set ip = If we assume that the back- 
ground evolution of the Universe is standard, with T~L 



and f n being exactly the same as in ACDM, Eq. (39) 
implies that the strength of the gravitational potential 
produced by density perturbations is modified by the 
factor Gn/G cos . When this modified potential is sub- 
stituted into the matter equations of motion (which have 
the standard form), it leads to a different growth rate of 
perturbations with respect to the ACDM case. For in- 
stance, a straightforward calculation shows that during 
the matter dominated epoch, the density contrast grows 
according to the modified power-law (cf. [12]) 



S oc a ^~ 1+ ^ 1+2AGN/Gcos) . 



(40) 



Notice that for small values of the parameters (a, /3, A) 
the anomalous growth is proportional to 



r, 1 = ^ + °( a ) , 



where we have defined 



a + P + 3A 



(41) 



(42) 



For E = 0, the effect of modified self- gravity is strongly 
suppressed. Hence we expect cosmological bounds on (a, 
/?, A) to be weaker along this degeneracy direction. In- 
cidentally, in the Einstein-aether case E is required to 
vanish (or, rather, be extremely small) by the PPN con- 
straints, see pbj). This explains why the cosmological 



bounds on Einstein- aether theory are rather mild [TTj fTB]. 



Model 


a 


P 


A 


E 


enhanced gravity 


0.2 





0.1 


0.5 


shear 


0.05 


0.25 


-0.1 






TABLE I: Parameters for the enhanced gravity and shear 
models. 



On the other hand, in the khronometric model the PPN 
constraints are compatible with E ^ and the influ- 
ence of Lorentz violation on cosmological perturbations 
is more pronounced. 

The second effect is understood from the tracefree part 
of the (ij) Einstein equation, 



k 2 ^ - <t>) 



UirG a 2 ^ (fi+Pifa + pk 2 ( X + ZHx)- 



(43) 



other 



Since the khronon field introduces a preferred direction 
(39) u \i m space-time, it may generate anisotropic stress. This 
effect (proportional to 0) is accounted for by the last term 
in the previous equation. Generally speaking, adding 
anisotropic stress amounts to increasing the viscosity of 
the cosmic fluid, and leads to a damping of small-scale 
perturbations. 

The third effect comes from gravitational interactions 
between ordinary matter species and the scalar fields x 
and £. This interaction, described by the second line in 
Eq. (37), may play an important role under the condi- 



tion that the fields cluster and form sufficiently dense 
clumps. This might be the case through a mechanism 
described in Ref. [4 . The mixing between the dark en- 
ergy perturbation £ and the khronon x ( see Eqs. (21)) 
gives rise to a mode whose sound speed vanishes in the 
limit of small momentum k. This property implies that 
the effective pressure associated with the mode is small, 
and that the density perturbations (Sp x , 5p^) can be am- 
plified efficiently by gravitational collapse. Hence, the 
OCDM model features clustering dark energy. A semi- 
quantitative analysis of this effect was performed in [4], 
showing that structure formation is enhanced at small 
comoving momenta (large wavelengths), k < H a . Unfor- 
tunately, in this range, the quality of cosmological data 
is rather poor, and this effect does not play a significant 
role in actual observational constraints on the model. 

To illustrate how the above effects impact observable 
quantities, we study numerically the evolution of cosmo- 
logical perturbations in two reference models using the 
Boltzmann code CLASS. We will call them the enhanced 
gravity and the shear models. The corresponding param- 
eter values are listed in Table [TJ Clearly, in the enhanced 
gravity model we keep the effect (i) while switching off 
the effect (ii); in the case of the shear model the situa- 
tion is opposite. The effect (Hi) is present in both models 
but we will see that it is always subdominant. Note that 
the values in Table [I] have been chosen very large in order 
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to make the modifications visible on the plots. However, 
these values are excluded by current data (see Sec. |V|). 



A. Effects on the CMB 

In this subsection, we describe the changes induced by 
Lorentz violation in the CMB temperature anisotropy 
spectrum. We consider 0CDM with the two reference 
sets of parameters listed above and compare the results 
to ACDM. To highlight the changes, all simulations are 
performed with adjusted initial conditions such that, 
in the limit k — » 0, the gravitational potential xp is the 
same for all three models. For the standard cosmological 
parameters, we choose the following values: n s = 1, 
h = 0.7, n b = 0.05, tt cdm = 0.25, A s = 2.3 x HT 9 , 

%reio — 10. 

Enhanced gravity model: In Fig. [I] (top panel) we can see 
that at the time of decoupling, the gravitational potential 
in this model is enhanced (in absolute value) compared 
to ACDM. This is due to an increase in the perturba- 
tion growth, governed by E. On the same panel one 
observes two changes in the solution for the photon tem- 
perature perturbations 7 = 5 1 /A: a shift of the peaks of 
oscillations towards higher momenta, and a shift of the 
zero point of oscillations. These effects can be qualita- 
tively understood from the combination of the modified 
Poisson equation (39) with the equations of motion of 



the photon-baryon plasma before decoupling. The latter 
have a standard form, and in the tight coupling approx- 
imation they reduce to a single master equation, 



e. 



R 



R^ 



-V- 



R 



R" 



(44) 



where R = |^ encodes the baryon-to-photon density ra- 
tio, and c s = (3(1 + i?)) -1 / 2 is the sound speed of density 
waves in the plasma in the absence of gravity. According 
to Eq. (|39|), the first term on the r.h.s. of (44) contains 



a contribution proportional to 7 . This contribution is 
positive and larger than in ACDM for E > 0. It ef- 
fectively decreases the speed of sound in the plasma at 
the moment when each mode enters inside the horizon 15 , 
which translates into a shift of the peaks, observed in the 
top panel of Fig. [I] Next, the zero point of the acoustic 
oscillations is given by the value of — ip/Sc 2 at decoupling. 
As already mentioned, \ip\ is larger in the enhanced grav- 
ity model, leading to a shift in this zero point further 
away from 7 = 0. Finally, the amplitude of the acous- 
tic oscillations around the zero point is slightly smaller 
in the enhanced gravity model than in ACDM. 



15 Strictly speaking, the Poisson equation ( |39| is valid only for sub- 
horizon modes. However, it is sufficient for our qualitative argu- 
ment. 



ACDM 

enhanced gravity Q y --- 
shear Q — 
ACDM -\|/ 
enhanced gravity -\\f - 

shear -\|/ — 




ACDM SW 
enhanced gravity SW 
shear SW 




k in h/Mpc 

FIG. 1: Top: photon temperature perturbation 7 = £ 7 /4 
and gravitational potential — ip at decoupling, for the three 
models under scrutiny. Bottom: Sachs- Wolfe contribution, 
given by B 7 + ip (i.e. by the difference between the two curves 
above) . 



The above features affect the Sachs- Wolfe (SW) contri- 
bution to temperature anisotropics. In the bottom panel 
of Fig. [Tj we plot this contribution, given by B 7 + ip at 
recombination. One clearly sees the shift of the peaks 
and notices that the even peaks are suppressed while the 
odd ones are almost constant (with the notable exception 
of the first peak, which is suppressed). This is due to the 
competition between the different effects described in the 
previous paragraph. 

In the anisotropics observed today, the SW contribu- 
tion is supplemented by those coming from the Doppler 
and Integrated Sachs-Wolfe (ISW) effects. The decom- 
position of the Ct spectrum in terms of these effects is 
presented in Fig. [2] for different models. The compari- 
son of the enhanced gravity model with ACDM is shown 
in the top-panel. In the SW contribution (as well as in 
the total spectrum), we observe the expected suppres- 
sion of the first, second and fourth peaks, while the third 
peak amplitude is roughly unchanged. On small angular 
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scales, the peaks are further suppressed by Silk damp- 
ing. Indeed, due to the shift in the phase of oscillations, 
they correspond to smaller physical scales at recombina- 
tion, that are more affected by diffusion damping. The 
Doppler effect, which depends on 6 7 at recombination, 
is also modified. Finally, a prominent feature clearly vis- 
ible on the plot is the significant enhancement of the 
ISW contribution in the range 10 < I < 100, i.e. be- 
tween the regions usually affected by the early ISW effect £ 
(100 < I < 200) and the late ISW effect (2 < I < 10). I 
The ISW effect is proportional to the time derivative of 
the gravitational potential. In the ACDM model, the po- 
tential varies only during the epochs of radiation and A 
domination. However, in the enhanced gravity model, the 
growth of density perturbations entails a slow increase of 
the gravitational potential also during the matter domi- 
nated era, enhancing the ISW effect on a wide range of 
scales. 

All in all, we conclude that the enhanced gravity model 
produces significant modifications in the spectrum of 
CMB anisotropics. The pattern of these modifications 
is quite specific, and apparently not degenerate with the 
effects of standard cosmological parameters. q 



Shear model: We recall that in this model, the \ field 
generates some anisotropic stress and contributes to the 
shear of the perturbed metric, as described by Eq. (43). 



The presence of shear tends to smooth out metric per- 
turbations on scales smaller than the sound horizon asso- 
ciated with the sound speed of the x field. Note that for 
parameters of the particular shear model studied here, 
the field x is superluminal 16 (c x = a/3) and the suppres- 
sion appears already on super-Hubble scales. 

This is indeed observed on the top panel of Fig. [TJ 
where the gravitational potential is clearly smaller (in 
absolute value) compared to ACDM. It also exhibits no- 
table wiggles caused by oscillations in the %-field [4 . The 
fact that c x is much larger than the photon-baryon sound 
speed explains the shift between the phase of the oscil- 
lations seen in ip and in 7 . The suppression of ip shifts 
the zero-point of the oscillations in the photon tempera- 
ture 6 7 . On the other hand, we do not see any shift in 
the positions of the peaks, which is compatible with the 
previous discussion: the self-gravity of radiation is not 
modified in this model. 

For fixed initial conditions 7 (to), the amplitude of 
the acoustic oscillations depends crucially on boosting 
effects, imprinted around the time of Hubble crossing, 
and caused by the three gravitational driving terms on 
the right-hand side of Eq. (44). In the shear model, the 
amplitude of acoustic oscillations is damped as a conse- 
quence of smaller metric fluctuations and reduced grav- 
itational boosting. This translates into an overall sup- 




16 As pointed above, this does not present any inconsistencies in 
theories without Lorentz invariance. 



FIG. 2: Temperature anisotropy spectrum (solid) and its 
decomposition in terms of Sachs- Wolfe (dotted), Doppler 
(dashed) and Integrated Sachs-Wolfe (dot-dashed) contribu- 
tions. For clarity, we do not show cross correlations between 
these contributions. Thick black lines represent the ACDM 
model, while thin blue lines are used for the two GCDM ref- 
erence models. 



pression of the SW effect visible in the bottom panel of 

Fig.H 

On Fig. [2| we see that the Doppler and ISW contri- 
butions to the total temperature spectrum Ci are also 
lower in the shear model than in ACDM. The net result 
is a uniform suppression of all peaks. One may expect 
that this effect could be compensated, at least partially, 
by a rescaling of the initial amplitude of perturbations. 
This suggests that pure shear models might be less con- 
strained than enhanced gravity ones. 

Our shear model is similar to the one studied in [IT] , 
where it was claimed that the dominant effect on the 
CMB comes through the ISW. Our analysis demonstrates 
that the changes in the SW and Doppler contributions 
are equally important for this model. 
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B. Effects on the matter power spectrum 



to the actual value, 



Another observable affected by the behaviour of cos- 
mological perturbations in 0CDM is the matter power 
spectrum. We restrict the discussion to low enough 
Fourier modes, k < 0.1 h/Mpc, for which perturbations 
are still very close to the linear regime. The study of 
non-linear clustering in this model is beyond the scope of 
the present paper. 

The definition of what one calls the "matter power 
spectrum" is not so obvious in 0CDM. This model con- 
tains additional components that contribute to pertur- 
bations of the total energy density. One may wonder 
whether they must be included in the calculation of the 
power spectrum. In general, the answer to this ques- 
tion is yes. Indeed, the existing observations fall into 
two categories. The first category probes directly metric 
perturbations (themselves related to total density fluctu- 
ations): this is the case e.g. for cosmic shear surveys or 
CMB lensing measurements. The second category mea- 
sures the clustering of compact objects like galaxies, ha- 
los or clusters. On large scales, these objects are known 
to trace linearly the underlying gravitational field [36] . 
Thus, it appears reasonable to define the matter power 
spectrum using the Poisson equation 17 



k 2 (j) = -47rG N a 2 5p t 



(45) 



Note that we have used here the locally determined value 
of the Newton constant introduced in Eq. (|8|. Comparing 
Eq. ( |45| ) to Eq. ( |37| ) and neglecting the anisotropic stress 
which is very small on relevant scales, we see that the 
contribution of ordinary matter to Sp tot coincides with 
the standard definition 5p ni while the contributions of 
the khronon and £ fields read 



S Px 



ak 2 



Sir a 2 G 
a 



87ra 2 G cog c|r a 



n(i - B) X ) 



(46) 



We stress that all quantities here are taken in the New- 
tonian gauge. 

To obtain the matter power spectrum, the density per- 
turbation Sptot must be divided by the total matter den- 
sity (dark matter plus baryons). The latter is determined 
from the background cosmological solution. In the case 
of standard gravity, one would again use the local New- 
ton constant Gat to infer the density from the geometry. 
Because in the OCDM model the gravitational constant 
in the Friedmann equation is different, the effective mat- 
ter density found in this way is renormalized compared 



Peff 



G c 



G 



N 



'-{pcdm + Pb) 



(47) 



Thus we arrive at the following formula for the power 
spectrum, 



p(k) 



G N V (\5ptot(k)f 



G C os J (pcdm + Pb) 



2 ' 



(48) 



Let us discuss the imprint of our two reference OCDM 
models on this spectrum. 

Enhanced gravity model: The left panels of Fig. [3] show 
the time dependence of various contributions to Sp to t in 
this model, for three values of the momentum k. As 
expected, inside the Hubble radius, matter density per- 
turbations grow at a larger rate and are enhanced com- 
pared to the ACDM case. The contribution of the x 
and £ fields, though small, exhibits some interesting fea- 
tures. On super-Hubble scales, the ^-density perturba- 
tion rapidly grows with time. This is due to the mixing 
between \ and £ discussed in Ref. [4 , which gives rise 
to a mode with vanishing sound speed at low momenta, 
and allows for clustering of the dark energy. In the case 
of modes crossing the Hubble scale around the current 
epoch, 5p£ even becomes comparable to the matter con- 
tributions at the present time. However, at shorter scales, 
the mixing between \ an d £ disappears and the speeds 
of sound of these components become non negligible 18 . 
This leads to damped oscillations of these fields, clearly 
visible in the lower left panel of Fig. [3] 

The ratio of the power spectra in the enhanced gravity 
model and in ACDM is presented in Fig. [4] The acceler- 
ated growth of matter density perturbations translates 
into a significant scale-dependent enhancement of the 
power spectrum on scales that are below the Hubble 
radius today (i.e. with k > 0.0003 fo/Mpc). The curve 
exhibits small wiggles due to a shift in the position of 
the peaks of baryon acoustic oscillations (cf. the shift 
of CMB peaks discussed in Sec. IV A). On sub-Hubble 



scales, the contribution of 5p x and 5p^ to the total 
matter power spectrum is negligible 19 . Finally, let us 
point out that the curve in Fig. [4] must be taken with 
a grain of salt at k > 0.1/i/Mpc, where non-linearities 
become important. 

Shear model: Our definition of the matter power spec- 
trum refers to sub-Hubble wavelengths only, and also uses 
the assumption that <j) and ip are equal in Eq. (37), which 



Strictly speaking, this definition makes sense only for modes well 
inside the horizon, where the Poisson equation is valid. However, 
this is sufficient for our study, because the power spectrum is 
actually measured for such modes only. 



We remind that in this work we consider the regime where c x , 
ce are of order one. If instead these velocities were very small, 
the fields x an d £ could in principle cluster also on scales much 
smaller than the Hubble radius. 

As mentioned above, this conclusion could change in a scenario 
where a tiny value of c x or c@ would be assumed. 
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FIG. 3: Evolution of the transfer function of various energy density perturbations 5pi (arbitrary units) in the Newton gauge, 
for the enhanced gravity model (left) and the shear model (right), and for the three wavenumbers k = 10 _7 Mpc _1 (top panels), 
fc = 6x 10 _4 Mpc _1 (middle panels) and k = O.lMpc -1 (bottom panels). For comparison, we show the evolution of 5pb and 
Spcdm in the ACDM case. We also indicate the conformal time of Hubble crossing for the two largest wavenumbers, and the 
conformal time of radiation to matter equality. 



leads to Eq. (45). Hence we need to check whether the 



anisotropic stress (responsible for a possible difference 
between the two metric perturbations, see Eq. (|43|)) van- 
ishes inside the Hubble radius in the shear model. It is 
well-known that the anisotropic stress of neutrinos and of 
decoupled photons decays inside the Hubble radius, due 
to free-streaming. This conclusion also applies to the 
anisotropic stress of the khronon, due to the dynamics 
of the x field. For instance, during matter domination, 



X oscillates with an envelope decaying like r -3 on sub- 
Hubble scale, implying that its averaged energy density 
decays like Sp x oc r -8 ex a -4 . The contribution of the 
khronon anisotropic stress to the difference (c/)—ip) is then 
oscillating with an envelope proportional to r -4 oc a~ 2 . 
At the same time, <p and ip are almost constant. Hence 
the relative difference between <fi and ip quickly becomes 
negligible, and our definition of the matter power spec- 
trum on sub-Hubble scale is applicable. 
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k (h/Mpc) 

FIG. 4: Ratios of the matter power spectra in the two refer- 
ence 9CDM models and in ACDM at redshift z = 0. 



If we normalize perturbations to the same initial value 
of the evolution of (5p c dm, Spb) is identical above the 
Hubble scale in the ACDM and shear models. However, 
in the shear case, the anisotropic stress of the khronon 
cannot be neglected soon after the time of Hubble cross- 
ing. This leads to a smoothing of metric perturbations 
and to a suppression of Sp c dm and 5pb that is clearly vis- 
ible in the lower right panel of Fig. [3j Well after Hubble 
crossing, the khronon anisotropic stress becomes negligi- 
ble and (Spcdrn, $pb) evolve like in ACDM, but with a 
constant offset coming from the suppression experienced 
soon after Hubble crossing. 

This suppression is seen better in Fig. |4j For the cho- 
sen parameter values, the effect of the shear model is 
smaller than for the enhanced gravity model. In order 
to see clearly both effects with the same scale, we mul- 
tiplied the difference between the power spectra in the 
shear model and in ACDM by 20. Inside the Hubble ra- 
dius, the suppression of the power spectrum in the shear 
model is almost scale-independent. When approaching 
the current value of the Hubble scale, the suppression is 
reduced by the counteracting effect of additional density 
fluctuations in the £ field. However, this reduction is only 
significant on scales that are too large to be observed with 
good precision (due to the sampling variance associated 
to a given survey). Besides, on such scales our definition 
of the matter power spectrum is no longer applicable 20 . 

We reach the same conclusion as in the CMB case: on 
observable scales, the effect of pure shear models can be 
mimicked by an overall reduction of the primordial fluc- 
tuation amplitude, which suggests that these models are 



Like in the enhanced gravity model, we note that in a scenario 
with tiny values of cq or c x , density fluctuations in the £ or % 
field could in principle be significant on much smaller scales. 



more weakly constrained than enhanced gravity models. 



V. COMPARISON WITH CURRENT DATA 

We will now compare the khronometric OCDM model 
to CMB and LSS data, using the parameter inference 
code Monte Python 21 [20]. For the CMB, we use here 
WMAP 7-year data [37], and SPT data from 2008 and 
2009 [38 . Updating our analysis with recent Planck data 
[39] could give a small improvement on parameter con- 
straints, without changing their order of magnitude. For 
LSS, we rely on galaxy power spectrum data from the 
WiggleZ redshift survey [40]. We choose to perform our 
runs in the synchronous gauge. 

We include in the fit eight free cosmological parame- 
ters, which are the usual six free parameters of the min- 
imal flat ACDM model, plus f3 and j3 + A (this combina- 
tion is chosen to facilitate the convergence of the chains) . 
The parameter a is fixed to 2/3, which is the condition 
for satisfying all bounds coming from Solar System tests 
in the khronometric model 22 . As far as cosmology is con- 
cerned, this condition is actually not necessary, but we 
implement it anyway in order to deal only with realistic 
models. We also vary three nuisance parameters describ- 
ing the foreground contamination of SPT data, and we 
marginalize over these parameters following strictly the 
approach of Ref . [38] . 

We show our results for the Bayesian minimum cred- 
ible interval of each parameter in Table [II] The one- 
dimensional and two-dimensional posterior parameter 
distributions are displayed in Fig. [5] (omitting the three 
SPT nuisance parameters for clarity). 

The OCDM parameters are found to be weakly cor- 
related with the six standard model parameters, and 
strongly correlated with each other. Isoprobability con- 
tours in the (/3, /3 + A) space have a complicated shape. 
They are very elongated along two directions of degen- 
eracy, corresponding to /3 + A ~ and /? ~ 3/2(/3 + A). 
Away from these two special directions, the contours are 
more regular, and closer to an ellipse centered on the ori- 
gin. The posterior probability peaks at (/?, /3+A) ~ (0, 0), 
showing that the data brings no evidence in favor of the 
OCDM model. The two directions of degeneracy can be 
interpreted as follows. 

First, the case j3 + A c± corresponds to E ~ for 
our choice of a. This is the region of pure shear models, 
discussed in the previous subsection. We already argued 
that the differences with ACDM are less important in this 



21 http://montepython.net 

22 We decided to focus only on the khronometric case, since from 
the comments after Eq. \A2\ it is clear that the bounds are 
stronger in this model than in the Einstein-aether case. This 
is confirmed by comparing our results with the constraints pre- 
sented in Ref. 
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TABLE II: Mean values and 68% confidence limits of the minimum credible interval of the parameters of the khronometric 
GCDM model. For f3 and /3 + A, we show instead 95% upper limits. 
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FIG. 5: Triangle plot showing the one-dimensional marginalized posterior distribution and the two-dimensional probability 
contours (at the 68, 95 and 99% confidence level) of the parameters of the khronometric GCDM model. 
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case 23 than in the enhanced gravity case corresponding to 
the orthogonal direction (3 = 0. This explains the factor 
four difference between constraints on j3 and on j3 + A in 
Table HI 

The other case has a less straightforward interpreta- 
tion. The relation j3 ~ 3/2 (/? + A) implies that the 
khronon has a squared sound speed c 2 ~ 1/3, similar to 
that of neutrinos and of tightly coupled photons. Hence, 
for these models, the khronon does not introduce a new 
characteristic scale, and the oscillations seen in the lower 
left panel of Fig. [3] may mimic those of photons or neutri- 
nos, thus lowering the sensitivity along this direction. We 
checked explicitly with the modified CLASS code that in 
this particular case, the sum of all effects described in the 
previous section (enhanced gravity, additional shear and 
extra clustering species) can be nearly cancelled at the 
level of CMB anisotropics by a variation of the standard 
cosmological parameters. This situation is analogous to 
the case of extra relativistic degrees of freedom, for which 
the impact of a small variation of N e ^ on the CMB can 
be partially compensated by a variation of other param- 
eters, leading to a well-known degeneracy between N e ^ 
and H gT]. 6CDM models with c 2 ~ 1/3 are not ex- 
actly equivalent to models with extra relativistic degrees 
of freedom, but their effect can be compensated in a sim- 
ilar way. This explains the weak correlation observed in 
Fig. [5] between standard parameters and 6CDM parame- 
ters. By comparing with a run based on CMB data only, 
we find that the inclusion of LSS data helps in break- 
ing this degeneracy (by limiting enhanced gravity effects 
on the matter power spectrum), but even in presence 
of WiggleZ data, the degeneracy appears very clearly in 
Fig. [5} 



VI. CONCLUSIONS 

In this work, we described the impact of the OCDM 
model on cosmological observables, on scales where linear 
cosmological perturbation theory is valid. This model is 
an alternative to the ACDM scenario in which the accel- 
eration of the universe expansion is caused by a dynam- 
ical scalar field O. An important difference between this 
model and most quintessence models is that the action of 
the field © is naturally protected from large ultraviolet 
corrections. In addition, the 6CDM model is embedded 
in a family of gravitation theories with broken Lorentz 
invar iance, which may have an ultraviolet completion in 



Under the condition a = 2(3, the velocity of the x-field is given 
by c£ = E/(6/3). Thus, the limit S = formally corresponds 
to vanishing sound speed of the x-component. One could expect 
that this would lead to an amplification of the x-perturbations 
and, as a consequence, to strong constraints on the model. How- 
ever, in practice, c x never becomes smaller than ~ 10 — 2 during 
the Monte Carlo sampling of the parameters and the clustering 
of x does not appear to be significant. 



the framework of Hof ava gravity. The latter theories are 
characterized by a time-like dynamical vector field u 11 
which defines a preferred time direction. Typical exam- 
ples are the Einstein-aether and khronometric theories. 
The OCDM model is the simplest extension of these the- 
ories to include a mechanism for cosmic acceleration. Un- 
der simple constraints on its parameters, it passes all lo- 
cal tests of gravity: Solar System dynamics, gravity wave 
emission, black hole structure. We have shown that the 
0CDM model may produce observable effects on cosmo- 
logical scales, constrained by current cosmological data, 
but still potentially detectable with future ones. 

The background evolution in OCDM is identical to 
ACDM. The only potential difference could come from 
a component playing the role of stiff matter, but in or- 
der to preserve primordial nucleosynthesis, we must as- 
sume that this component is too small for playing any 
role in late-time cosmology, and for affecting CMB and 
large scale structure observables. However, the evolution 
of cosmological perturbations is generically very different 
in the OCDM and ACDM models. We studied the evo- 
lution equations of scalar perturbations in OCDM and 
identified three different effects: (i) a renormalization of 
the matter contribution to the Poisson equation (i.e. a 
different self-gravity of matter perturbations), (ii) a new 
contribution to the anisotropic stress, and (Hi) the pres- 
ence of additional clustering degrees of freedom. The first 
two effects are generic for Lorentz violating gravitation 
theories based on a unit time-like vector, while the third 
one is specific to the dynamical realization of dark energy 
in OCDM. 

We implemented the equations of OCDM in the Boltz- 
mann code CLASS [T9] . in order to compute accurately 
the impact of the three effects (i) , (ii) , (Hi). We found 
that they affect the power spectrum of Cosmic Microwave 
Background anisotropics in a very particular way (shift 
in the position and amplitude of the peaks, and enhanced 
ISW), as can be seen in Fig. [2] Furthermore, they influ- 
ence the shape of the matter power spectrum (different 
amplitude and slope on observable scales, shift in the po- 
sition of baryon acoustic oscillations). These effects are 
shown in Figures [3] and [4) 

To derive constraints on the free parameters of the 
khronometric OCDM model (considering only parameter 
combinations satisfying Solar System tests), we used data 
from WMAP (7-year), from SPT (2008 and 2009 obser- 
vations), and from the WiggleZ redshift survey. We ran 
the parameter inference code Monte Python [20] and 
found the bounds displayed in Table [IT] and illustrated in 
Fig. [5j Quite remarkably, these bounds constrain devi- 
ations from general relativity to better than the percent 
level. They are stronger than those from several tests of 
gravitation, such as radiation damping of binary systems, 
or BBN. They are also stronger for the khronometric case 
studied in this paper than for the Einstein-aether model 
considered in Ref. [11 , due to the absence of the effect (i) 
in the latter case. The vanishing of (i) in the Einstein- 
aether theory (once Solar System bounds are imposed) 
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seems to be simply a coincidence. 

The fact that our bounds are dominated by the ef- 
fect (i) - modified self-gravity of the matter perturba- 
tions - suggests that they may have a wider application. 
The bound on the combination (3 + A in Table II can be 
cast into the constraint on the discrepancy between the 
gravitational constants appearing in the Newton law and 
in the Friedmann equation, \Gn/G cos — 1| < 0.018 at 
95% confidence level. Though in our model we cannot 
completely disentangle (i) from other effects, we believe 
that the above constraint will apply, at least by order of 
magnitude, to any theory predicting a time-independent 
discrepancy between Gn and G cos . It represents an im- 
provement compared to the bounds existing in the liter- 
ature nana]. 

Our analysis can be extended in several ways. First, 
we have only used the scalar sector of the theory. For the 
Einstein- aether case, the theory contains also propagat- 
ing vector modes, whose influence on B-type CMB polar- 
ization could provide further observational tests [I3j [15] . 
Second, the analysis can be updated with CMB data from 
Planck, and in the future, with large scale structure data 
from DES, LSST or Euclid. It would also be extremely 
interesting to go beyond the linear regime and try to 
understand the consequences of OCDM for non-linear 
structure formation. Finally, to completely characterise 
possible deviations from Lorentz invariance in the con- 
text of cosmology, one can consider the option that this 
symmetry is violated also in the dark matter sector [27] . 
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the perturbations is provided by the growing adiabatic 
mode. In this Appendix we determine this mode ana- 
lytically on super-Hubble scales using the equations of 
Sec. UTTBl 

Taking into account that a ex r during radiation dom- 
ination (r stands for conformal time), and keeping only 
leading order terms in an expansion in (k/H) ~ (fcr), we 
can simplify Eqs. (21a) and (21b): 



x- 

I- 



2B 



r 



^^+^+^7 = 0, (Al) 

^cos >ot ^ ™ 

(A2) 



On the other hand, Eqs. (24), (28) can be rewritten as 
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2k 2 G cos t 2k 2 r 2 G c 
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(R V 5 V + (1 - R„)6 7 ) 

a Go Ha L (A3) 



h 6 

h+- + (RvSv + (1 - Rv)S 7 ) 

T T 7 



G 



G V r 



(A4) 



where R v is the ratio of the density of neutrinos to 
the total radiation density. Since we have not modi- 
fied the matter sector of the theory, the solutions for 
matter perturbations in the super-Hubble regime coin- 
cide with Eqs. (92) of [31]. Note that we have left aside 
the remaining pair of Einstein equations ([26]) ? (30), that 



are redundant and automatically satisfied by the solution 
dA7b— dA10). 



At initial time To, the 0-field contribution to the total 
energy density is negligible. It is natural to expect that 
its fluctuations can also be neglected, which c orres p onds 
to dropping off the terms containing £ in Eqs. ( Al ), (LA3I) 



and (A4). Let us find the precise conditions under which 



this is possible. The non-decaying solution of (A2) is 



(A5) 



where £o is a cons tant i mply ing that the £ con tribution 
drops out of Eqs. ([A3]) , (lA4|). As for Eq. ([All), we find 



that the £-term is negligible provided the following con- 
dition is satisfied: 



£o < -— h 

-tlrv 



(A6) 



Appendix A: Initial conditions 

The initial conditions for the system of differential 
equations describing the evolution of a given Fourier 
mode k are set in the radiation era, at a moment To 
when the wavelength associated to k is well outside the 
Hubble scale, i.e. when k <C T~L. As discussed in [4j[12], 
under broad assumptions, the dominant contribution to 



where we have used the expression (|A7a|) for h and have 
assumed ce,c^ 



1. This inequa 



ity is the condition 



for having negligible isocurvature perturbations in the £- 
field. It is likely to be satisfied under plausible assump- 
tions about the origin of the primordial ^-perturbations 
[4]. However, this statement must be taken with a grain 
of salt: a detailed theory of the O-field dynamics starting 
from the inflationary epoch is required to put it on the 
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solid ground. We leave developing such a theory for the 
future. In the present work we will set £o to zero in the 
numerical simulations. 

For other fields we use the ansatz 24 (cf. [31 J: 



h = h r + hir , r\ = T] + 771 r , 



(A7a) 



tf 7 = £ 7 qt 2 + ^ 7 ir 4 , b v = S u0 r 2 + ^it 4 , (A7b) 



6/ 7 ot 

2 



or 

3 



(A7c) 
(A7d) 



From the matter equations (Eqs. (92) of [31]) we obtain 
the standard relations: 
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To fi n d th e rela tions between the fields /io, ^1,^1 we use 
(A3), (A4) and (Al). In the leading order the first equa- 



tions is satisfied identically, while the second yields 

2G 7 



k 2 G c 



(A9) 



Considering the subleading order in (A3), (A4) and the 
leading order in (Al) we obtain 
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It is now straightforward to check that the obtained so- 
lution satisfies the remaining Einstein's equation (26), 
(|30). 



To summarize, the adiabatic mode is given by 



Eqs. ( |A7| ) with the coefficients related by ( |A8| ), ( |A9| ), 
(A10) to the single constant ho setting the overall nor- 
malization. To these expressions we add the standard 
adiabatic initial conditions for baryons and dark matter 



cdm 



, e h 



(AH) 



-h 

— ^ 7 ' u cdm — VJ ; [y b — ^7 

We cross-checked the validity of our initial conditions 
using the numerical code. When fixing initial conditions 
at ro, we find that over an extended range of time above 
ro, the numerical solutions for (fo, 77, \) an d for all density 
perturbations remain equal to their analytic expressions 
in Eqs. (A7)-(A10). Hence, the above initial conditions 



correctly describe the attractor solution of the full sys- 
tem of eq uation s that is compatible with the adiabatic 
condition (All) for matter fields. At leading order, the 



attractor solution for the field £ is given by 



(cokr) 2 



42H„ 



(A12) 



and evolves proportionally to r 6 . Since this solution is 
too small to back-react on other fields, it makes no dif- 
ference to fix the initial £ to the above solution or to zero 
(in the latter case, £ reaches very quickly the attractor 
solution) . 

When deriving initial conditions, we assumed full ra- 
diation domination, and neglected all corrections related 
to the contribution of non-relativistic matter to the ex- 
pansion rate. Hence our initial conditions are accurate 
only if imposed at very early time. We decreased ro in 
the code until getting very stable results. This is the 
case for the value ro = 10 _2 Mpc, that we adopted in all 
simulations. 
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